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ABSTRACT 

Considered  methods  of  solving  Kolmogorov  -  Fisher  type  task  of  reacting  with  diffusion,  in  the  two-dimensional 
case,  obtained  invariant  properties  of  solutions  and  two  way  estimation  of  the  solution.  Given  are  the  results  of  numerical 
experiments  for  various  values  included  in  the  equation  parameters,  in  the  two-dimensional  case. 
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1.  INTRODUCTION 

Currently,  studies  are  on  various  properties  of  solutions  of  nonlinear  problems  in  large  number  of  works,  and 
identifying  all  new  and  newer  properties  of  their  solutions.  Solution  of  nonlinear  boundary  value  problems  are  always 
accompanied  by  considerable  difficulties,  because,  to  solve  them  in  an  analytical  form  is  possible,  only  in  exceptional 
cases  and  for  establishing  new  properties  of  the  solutions  require  painstaking  research.  Standard  methods  of  research  on 
nonlinear  problems,  depending  on  the  type  of  nonlinearities  are  not  yet  available.  Therefore,  in  each  case,  to  study  the 
properties  of  solutions  have  resorted  to  various  exact  and  approximate  methods,  in  leading  universities  and  research  centers 
around  the  world,  including  Oxford  University,  University  of  Cambridge,  Department  of  Mathematical  Sciences  of 
University  of  Liverpool,  Department  of  Applied  Mathematics  of  the  University  of  Leeds  Faculty  of  Biological  Sciences  of 
the  University  of  Leeds,  Department  of  Engineering  Mathematics  and  School  of  Biological  Sciences  of  University  of 
Bristol  (UK),  Department  of  Mathematics  of  Stanford  University,  Consortium  of  the  Americas  for  Interdisciplinary 
Science  and  Department  of  Physics  and  Astronomy  of  University  of  New  Mexico,  Department  of  Biology  of  University  of 
Louisiana,  Departments  of  Entomology  and  Biology,  Pennsylvania  State  University  (USA),  Ecole  des  Flautes  Etudes  en 
Sciences  Sociales,  Institut  de  Mathematiques  de  Toulouse,  Universite  Paul  Sabatier  (France),  Universit'a  degli  Studi  di 
Padova  Dipartimento  di  Matematica  (Italy),  Departamento  de  Flsica,  Facultad  de  Ciencias  Fisicas  y  Matematicas, 
Universidad  de  Chile  Complex  Systems  Group,  Facultad  de  Ingenierfa  y  Ciencias  Aplicadas,  Universidad  de  los  Andes 
(Chile),  Centro  Ato  mico  Bariloche,  Instituto  Balseiro  and  CONICET  (Argentina),  Department  of  Theoretical  Ecology, 
Biology  Centre  all  ASCR,  Institute  of  Entomology  (Czech  Republic),  School  of  Sciences,  Jimei  University,  Xiamen 
(China),  Moscow  State  University,  Institute  of  applied  mathematics,  Institute  of  theoretical  and  experimental  Biophysics, 
Tomsk  State  University  (Russia),  national  University  of  Uzbekistan,  Samarkand  state  University  (Uzbekistan). 

Review  of  international  research,  based  on  the  works  of  Oxford  University,  University  of  Cambridge,  Department 
of  Mathematical  Sciences  of  University  of  Liverpool,  Department  of  Engineering,  Mathematics  and  School  of  Biological 
Sciences  of  the  University  of  Bristol,  Department  of  Mathematics  of  Stanford  University,  Consortium  of  the  Americas  for 
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Interdisciplinary  Science  and  Department  of  Physics  and  Astronomy  of  the  University  of  New  Mexico  [1-7]. 

Review  of  foreign  scientific  research  shows  that,  since  the  work  of  Turing,  the  mathematical  model  of  the 
formation  and  propagation  of  nonlinear  waves  and  processes  of  a  structured  self  -  organization  in  physical,  chemical, 
biological,  and  social  systems  such  as  reaction-diffusion,  where  members  describe  nonlinear  kinetics,  and  transport 
processes,  represented  by  isotropic  diffusion  (Oxford  University,  University  of  Cambridge).  However,  in  many  systems,  no 
less  important,  more  complex  mechanisms,  diffusion  is  a  nonlinear,  anisotropic  and  cross-diffusion.  Most  studies  have 
examined  direct  diffusion  (self-diffusion)  described  by  the  equations,  the  diffusion  coefficients,  which  are  constant  values 
(Ecole  des  Hautes  Etudes  en  Sciences  Sociales,  Institut  de  Mathematiques  de  Toulouse,  Universite  Paul  Sabatier  (France), 
Universit'a  degli  Studi  di  Padova  Dipartimento  di  Matematica  (Italy)).  In  this  class  of  systems,  the  processes  of  formation 
of  spatial-temporal  structures  are  determined  by  the  diffusion  coefficients  and  the  specific  forms  of  the  kinetic  function  of 
the  reaction  process.  The  case  when  the  diffusion  coefficients  are  not  constants,  but  depend  on  dynamic  variables, 
corresponds  to  the  nonlinear  diffusion.  Examples  of  nonlinear  diffusion,  in  the  processes  of  mass  transfer  in  porous 
environment,  as  well  as  in  population  models  (Consortium  of  the  Americas  for  Interdisciplinary  Science  and  Department 
of  Physics  and  Astronomy  of  the  University  of  New  Mexico).  Mathematical  models  with  diffusion  coefficient  depending 
on  the  density  of  bacteria,  describes  the  processes  of  formation  of  complex  spatial  structures  with  the  growth  of  bacterial 
colonies.  Regimes  with  aggravation  in  the  spatial  localization,  in  open  dissipative  systems,  described  by  models  with 
nonlinear  diffusion  (Institute  of  applied  mathematics,  Institute  of  theoretical  and  experimental  Biophysics,  Tomsk  State 
University). 

2.  STATEMENT  OF  THE  TASK 


Consider  a  two  dimensional  problem  of  Kolmogorov  -  Fisher  type  for  the  reaction  diffusion 
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In  the  domain  D  =  £2x(0,T),  O  cz  R2 ,  0  =  {—  ba  <xa  <ba,  <2  =  1,2}  with  the  initial  and  boundary 
conditions 


u(0,x )  =  u0  (x)  >  0, 


(2) 


U  —  jU(  X,  t )  ,  t  £  (0,  T)  ,  T  -  boundary  of  12  . 


To  solve  this  problem  (1),  (2)  uses  the  initial  approximation 
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In  this  case,  obtained  the  following  condition  for  the  localization  of  solutions  of  problem  (1),  (2):  T(°°)  <  +°°  , 
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Note  that,  the  way  to  build  initial  approximation  based  on  the  splitting  of  the  equation  (1)  [1]. 

This  task  in  the  case  m-0  h  k(t,x)-const  was  studied  in  several  papers  [2,3-5]. 

In  [6]  derived  the  invariant  properties  of  solutions.  As  in  [2]  for  the  models  of  Malthus,  Ferhulst,  Olli  obtained 
properties  of  the  finite  speed  of  the  perturbations.  In  [7]  established  properties  of  flash  localization. 

PROPERTIES  OF  INVARIANCE  OF  THE  SOLUTION 


Let's  show  that  the  function 

u(t,  x)  =  y/(t)w(z{t),  x),  (3) 


where  l//(t )  solution  to  the  equation  without  diffusive  part  of  equation  (1) 
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again  satisfies  the  equation  (1). 

In  fact,  after  putting  (3)  into  (1)  subject  to  (3),  (4)  it  is  easy  to  calculate  that  for  \v{t ,  x)  have  the  equation 
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where  w(r(f),  x)  new  unknown  functions  T(t)  -  to  be  determined  function. 
\//(t )  we  find  from  the  equation: 
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Substituting  (3)  into  (1)  and  choosing  T(t)  from - =  \ff  have: 
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where  Xj/x (t)=  k(t)y/^  °  ■ 

t 

However,  it  is  clear  that  due  to  (4)  lim  lj/(t)  =  1,  if  k\t)dt  exist.  Therefore,  we  can  assume  that  for  sufficiently 

/— >oo  J 

0 

large  t ,  \f/x  (?)  ~  k(t) .  ie  will  again  get  the  equation  (1).  Because  of  this,  we  call  the  function  y/(t)w(r(t),  x)  where 
y/(t)-  the  solution  of  equation  (4),  and  Vf(t(T ),  x)  solution  of  equation  (5)  invariant  of  the  equation  (1). 

For  the  solution  of  the  problem  (1),  (2)  fair 

Theorem  1 


LetO  <  u0  (x)  <  1,  ie  R,  £  =  |jc| / f1/2 . 

Then  for  the  solution  of  the  task  (1),  (2)  in  Q  —  {(/,  x)\  t  >  0,  X  6  R  '  }  takes  place 


two-way  estimation 


Proof 


1  j  a(t])d!j  1 

y/(t)(T  +  T(t)YNK2+aN) (a-^2)°<  u(t,  x)  <e°  (T  +  T(t)TN K2+aN) (a 

where  y/(t )  defined  above  function,  N-2. 


For  the  proof  of  the  theorem  first  get  the  estimate  from  above.  To  this  end,  in  (1)  will  replace 
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Therefore,  according  to  comparison  theorem  of  solutions  [1]  has  an  upper  bound 
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In  order  to  obtain  a  lower  bound,  we  apply  the  method  of  nonlinear  splitting  [1],  According  to  this  method  the 
bottom  the  solution  is  sought  in  the  form 

u(t,  X )  =  i/f(t)w_  (r(t),  x),  (9) 


Where  l//(t )  above-defined  by  formula  (4)  function. 
Then,  from  (5)  we  have 
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k{t)yP  n W_  (l  —  W_P  )  >  0,  in  Q,  if  constant  T  ~>  1 . 


Then  applying  the  comparison  theorem  of  solutions  [1]  because  of  (10)  have 

U(t,  x)  >  r(t)(T + T(t)y",atM\a 


with  (9)  proves  the  validity  of  theorem  1 . 

It  is  an  important  generalization  of  this  task  on  different  occasions.  In  particular,  it  has  not  been  studied  in  the  case 
of  heterogeneous  medium  (the  diffusion  coefficient  is  a  function  of  the  spatial  variable),  the  reaction  coefficient  depends 
on  time  or  has  a  more  complex  character. 
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=  jU(  X,  t )  ,  t  €  (0,  T)  ,  T  boundary  of  . 


To  solve  this  task  (11),  (12)  uses  the  initial  approximation 
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in  case 
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In  this  case,  obtained  the  following  condition  for  the  localization  of  solutions  of  task  (11),  (12):  T(°°)  <  +°°  , 

jk(r))dri 

q\  e°  drj  <  +°o  ■ 


Note  that,  the  way  to  build  initial  approximation  based  on  the  splitting  of  equation  (11)  after  bringing  it  to  the 
radially  symmetric  form  [  1  ] . 

Show  that  the  function 

u(t,  X )  =  if/(t)w{z(t),  x),  (13) 


where  y/{t)  solution  to  the  equation  without  diffusive  part  of  equation  (11) 

dy/ 
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=  k(t)lff(\-yp)  ,T.e.  y/(t)  = 


k(t)dt 

1  +  e  0 


(14) 


Again,  it  satisfies  the  equation  (11). 

In  fact,  after  putting  (13)  into  (11)  with  (13),  ( 14)  it  is  easy  to  calculate  that  for  w{t,  x)  have  the  equation 
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y/(t)  we  find  from  the  equation: 


(^-  =  kit)yri\-yrp),  Vit)  = 

dt 


i 

1 


V 


J 


. .  (It  g 

Substituting  (13)  into  (11)  and  choosing  Tit)  from  - =  \j/  we  get: 
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where  —  k(t)y/^  a  . 


However,  it  is  clear  that  due  (14)  lim  —  1 ,  if  j" k(t)dt  exist.  Therefore,  we  can  assume  that,  for  sufficiently 


o 

large  t,  y/\(t)  ~  k[t).  i.e  again  obtain  equation  (11).  Because  of  this,  we  call  the  function  y/(t)\\{T(t),  X )  where 
solution  of  equation  (14),  and  w(T(t),  x)  solution  of  equation  (15)  invariant  of  the  equation  (11). 


Below  we  study  the  effect  of  heterogeneity  in  m--2(N-l),  N  >  2,  where  N  dimension  of  the  space  in  the 
following  cases:  jB  <  (J  +  1 ,  jB  —  C 7  +  1  and  jB  >  (J  +  1 . 

Simulated  picture  of  the  evolution  process  is  on  the  computer  by  using  Matlab. 

4.  NUMERICAL  EXPERIMENT 

The  following  are  the  results  of  numerical  experiments  for  various  values  included  in  the  equation  parameters  in 
the  two-dimensional  case: 


te  [0,1];  x1  g  [-6,6];  x2  e  [-6,6];  n  =  1.1;  nXi  =  30;  nXi  =30;  nt  =100. 

In  all  considered  cases  the  proposed  approach  number  of  iterations  on  average  did  not  exceed  three  at  a  given 
precision  eps. 

Created  on  input  language  Matlab  the  program  allows  you  to  visually  trace  the  evolution  process  for  different 
values  of  the  parameters  and  the  data  (Fig.  [1-4]). 

File  Edit  View  Insert  Tools  Window  Help 

DcSH#  IA/  /  &  J3  O 


Figure  1:  Results  of  Numerical  Experiments  at:  X{  =  1;  X0  =  1 ;  O  =  1 ;  m  =  0.77 


Figure  2:  Results  of  Numerical  Experiments  at:  X1  =  2;  X2  =  2 ;  (J  =  2  ;  m  =  0.87 
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Figure  3:  Results  of  Numerical  Experiments  at:  X1  =  1;  X2  =  1 ;  <7  =  1;  m  =  1.77 


Figure4:  Results  of  Numerical  Experiments  at:  Xj  =  2;  X2  =2;  <J  =  2  ;m  =  1.87 
5.  CONCLUSIONS 

The  data  of  numerous  experiments  show  that  the  application  of  computer  technology  to  the  problems  of 
populations  for  a  wide  range  of  changes  in  input  parameters  allows  to  obtain  results  with  a  high  degree  of  accuracy  and 
with  little  expenditure  of  computer  time. 
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